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Abstract: S(pacelike)D-branes are objects arising naturally in string theory when Dirich- 
let boundary conditions are imposed on the time direction. SD-brane physics is inherently 
time-dependent. Previous investigations of gravity fields of SD-branes have yielded unde- 
sirable naked spacelike singularities. We set up the problem of coupling the most relevant 
open-string tachyonic mode to massless closed-string modes in the bulk, with backreaction 
and Ramond-Ramond fields included. We find solutions numerically in a self-consistent ap- 
proximation; our solutions are naturally asymptotically flat and time-reversal asymmetric. 
We find completely nonsingular evolution; in particular, the dilaton and curvature are well- 
behaved for all time. The essential mechanism for spacetime singularity resolution is the 
inclusion of full backreaction between the bulk fields and the rolling tachyon. Our analysis 
is not the final word on the story, because we have to make some significant approximations, 
most notably homogeneity of the tachyon on the unstable branes. Nonetheless, we provide 
significant progress in plugging a gaping hole in prior understanding of the gravity fields of 
SD-branes. 
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1. Introduction 



SD-branes in string theory were first studied by Gutperle and Strominger in ref. [Q]. They 
were introduced as objects arising when Dirichlet boundary conditions on open strings are put 
on the time coordinate, as well as on spatial coordinates. SD-branes are not supersymmetric 
objects, which makes them hard to handle but potentially very interesting. The boundary 
conditions for SD-branes imply that they live for only an "instant" of time, and so the 
worldvolume is purely spatial. SD-branes should not be confused with instantons, because 
they live out their lives in Lorentzian signature context. Recent discussion of the relation 



between SD-branes and instantons may be found in section 6 of ref. |1C]. 

SD-branes are especially interesting objects to study in the context of tachyon conden- 
sation, which will be the arena of our investigation. SD-branes are indeed inherently related 
to the general study of time dependence in string theory. One of the original goals of ref. [||] 
was in fact to seek examples of gauge/gravity dualities where a time direction on the gravity 
side is holographically reconstructed by a Euclidean field theory.^ 

There have been several investigations of SD-branes since they were introduced ||^, ^, ^, 
^, ^, |9|, Most of them involve taking the limit gg — > 0, the regime in which perturbative 
string computations can be done. As the tachyon rolls down its potential hill, there is a 



divergence in production of higher mass open string modes |q, 11, M. This divergence occurs 



before the tachyon gets to the bottom of the potential well, as it must because there are no 



perturbative open string excitations around the true minimum of the tachyon potential |12]. 
Also, the time taken to convert the energy of the rolling tachyon into these massive open 
string modes is of order^ ^{9s) ■ This analysis was done for a single SD-brane using CFT 



methods; analysis for production of massive closed string modes was also done |13, 14 1. 

One aspect of SD-brane physics has become clear: that the decoupling limit applied to 
SD-branes is not a smooth limit like it is for regular D-branes. In particular, as g'^ ^ the 
brane tachyon becomes decoupled from bulk modes, which were however the most natural 
modes into which the initial energy of the unstable brane should decay. Then, the endpoint of 
the rolling tachyon must include a somewhat mysterious substance called "tachyon matter" 
[ |l^ ] . Consideration of the full problem with gg finite would presumably eliminate the need for 
mysterious tachyon matter; this was in fact part of our motivation for this work. Regarding 
production of closed string massive modes, at very small gg, it seemed that there was some 
debate [|^, Q about the form of a divergence. The results from ref. |14] make clear that the 
divergence depends on the number of dimensions transverse to the decaying unstable brane: 
for unstable Dp-branes with p < 2 there was a need to invoke a cutoff to get a finite result. In 
any case, unstable brane decay should presumably be a physically smooth process for {gs,is} 
finite. 



^An attempt to find a realization for such a duafity is the dS/CFT correspondence |^ - see also ref. Q for 
an extensive list of references. 

^In later sections we will see that our calculation, which includes gravity backreaction, agrees with this in 
the sense that the time it takes for the tachyon to decay only slightly depends on the particular value of g^. 
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The point of view that we will be taking is to consider a system of N SD-branes, with 
QsN large, study the overall centre-of-mass tachyon, and couple it to bulk massless closed 
string modes. Obviously, it would be nice to understand the full problem including coupling 
to all massive open and closed string modes, but this is a hard problem beyond our reach. 
We will make a beginning here with a quantitative analysis involving only the lowest modes 
in each of the open and closed string sectors. We believe that our approach, while "lowbrow" 
by comparison to SFT computations, already shows some very interesting physics. 

The punchline of our paper will be this: we find nonsingular solutions for the evolution 
of the open string tachyon coupled to bulk supergravity modes. This plugs a gaping hole in 
our previous understanding [Q, ^, ^ of the supergravity fields arising from a large number 
of SD-branes. All previous attempts at describing SD-branes in the context of supergravity 
had found that the corresponding solutions were plagued with naked spacelike singularities. 
We find that resolution of these singularities is achieved in a conceptually simple way: by 
including full backreaction on the rolling tachyon. 

Our investigation can also be considered to shed light on the question of tachyon cos- 
mology [ p!^ ] including Ramond-Ramond fields, the effect of which was ignored in previous 
investigations. Tachyon cosmology itself may not yet provide realistic models for inflation, 
nonetheless, see recent work including, e.g., refs. [0, [l^. One of the reasons is that, in the 
low-energy actions used to describe tachyon cosmology dynamics, there is only one length 
scale — the string length. It would certainly be interesting if a mechanism generating a 
lower scale for inflation were found within this context. Also, the behavior of inhomogeneities 
during the later stage of the roll of the tachyon may be a general problem [^]. Tachyon 
cosmology involving brane-antibrane annihilation may be relevant only to a pre-inflationary 
period, but it is interesting to analyze the dynamics from the "top-down" perspective in string 
theory. 

The plan of our paper is as follows. We begin by reviewing in section ^ the previous work 
on gravity fields of SD-branes, and commenting on the nature of singularities found in the 
past. In section 2A we discuss pathologies of the solutions found in refs. @, |5|, where only 



supergravity fields were considered. In section we move to discussing ref. ||7| , in which an 
unstable brane probe was coupled to a d = 4 SDO-brane supergravity background; we general- 
ized their arguments but still find generic singularities in the probe approximation (exceptions 
are considered in Appendix In section ^ we set up the full backreacted problem of interest. 
The equations are naturally highly nonlinear, and since backreaction is essential we have no 
desire to ignore it or treat it perturbatively. We need to use a particular homogeneous ansatz 
to facilitate solution of the equations of motion, and we discuss implications of the ansatz. 
In section^ we demonstrate numerical solutions of our backreaction-inclusive equations, and 
interpret qualitative features found. In particular, we follow carefully the evolution of both 
the dilaton and curvature invariants. Lastly, in the Discussion section we summarize our 
conclusions, open issues, and directions of future work. 
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2. Singular supergravity SD-branes: Review 

In the paper Q, a small number of supergravity solutions, thought to be appropriate for 
a large number A'^ of SD-branes, were presented. Subsequently, the class of supergravity 
solutions was widened considerably, simultaneously by two groups (refs. and p). A later 
paper showed that these two sets of solutions were equivalent by matching boundary 
conditions on asymptotic fields and finding the coordinate transformation explicitly. 

2.1 Sourceless SD-brane supergravity solutions 

The convention of ref. which we will use, is that SDp-branes have (p+l) worldvolume 
coordinates. We call these y. There is also the time coordinate t, and the (8— p) overall 
transverse coordinates x. 

The most general SDp-brane supergravity solutions of refs. and |^ can then be written 
in string frame as follows. 



+F{ty 



-1/2 



.i=2 



Vow; w'^Hiji 



g(t)\-g-"*' 



m^r^^) , (2,1) 



C(P+i) = sin e cos 9 ^ dy^ A ■ ■ ■ A dy 
F{t) ^ ^ 



where 



and 

+(^)'-^ m^l-{j) (2.3) 

The supergravity equations of motion will be satisfied when the exponents satisfy the con- 
straints, 



P+l 



2 = 1 

p+l J _ 

~k + Y.h-^{H-G) = Q, (2.4) 

i=l 

4 

H + G-- = . 

7-p 

The general metric above is not isotropic in the worldvolume directions y. However, 
from a microscopic point of view, one expects that the supergravity solution would have an 
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isotropic worldvolume. Isotropy in the worldvolume will be restored in the above solution 
for the choice — ^2 = " ■ ■ = ~kp+i = ki + k = n. Unfortunately, the isotropy requirement 
excludes SD-brane solutions with regular Cauchy horizon. In fact, the curvature invariants 
associated to all isotropic solutions diverge at t = ibw and t = 0. Consequently, although they 
possess the right symmetries and "charge", these solutions do not appear to be well-defined. 

Besides, SD-branes should be represented by solutions with, roughly, three distinct re- 
gions: the infinite past with incoming radiation only in the form of massless closed strings, 
an intermediate region with both open and closed strings, and finally the infinite future with 
dissipating outgoing radiation in the form of massless closed strings. But the supergravity 
solutions of refs. ^, |5| cannot represent this process, because there are no rules for deciding 
how to go through the singular regions (see, however. Appendix ^ where we show that some 
anisotropic solutions avoid the pathology). 

Nevertheless, the isotropic solutions have a positive feature worth noting: they have the 
correct asymptotics at large time. In the limit that the functions a{t) and P{t) become trivial, 
part of the metric is simply the Milne universe: flat Minkowski space foliated by hyperbolic 
sections, 

p+2 

hm dsl^ = -dt' + t^dHl_p{x) + y2{dy'f. (2.5) 

i=l 

On the other hand, there are (at least) two reasons to suspect that the above solutions are 
not the final word in the SD-brane supergravity story. The first is that there are one too 
many parameters in the solution, by comparison to expectations from microscopies of SD- 
branes [^. A possible explanation for this may be as follows. In the rest of our paper, we 
will be showing that the full coupling/backreaction between the open-string tachyon and the 
closed-string bulk modes is crucial for resolution of spacetime singularities. It is possible that 
the freedom in the supergravity solutions may correspond to a freedom in picking boundary 
conditions for the rolling tachyon — the coupling to which was not included in refs. |]l|, ^, ^. 

The most noticeable negative feature of the above solutions is that the isotropic solutions 
are nakedly singular. Quite generally, nakedly singular spacetimes arising in low-energy string 
theory come under immediate suspicion, even though they are solutions to the supergravity 
field equations. No-hair theorems are usually what we rely on in order to be sure that we 
have the unique supergravity solution, but no-hair theorems are never valid for solutions with 
naked spacetime singularities. It is worth noting that it has been shown with an explicit coun- 
terexample that even no-hair theorems themselves fail for black holes in d = 5 with mass 
and angular momentum — and hence the idea of no-hair theorems in all higher dimensions 
is under suspicion. (Nonetheless, with particular assumptions about field couplings, no-hair 
theorems can be proven for static asymptotically flat dilaton black holes |2^. Also, unique- 
ness of the swper symmetric rotating BMPV black hole in d = 5 has been proven |22|.) 
Even if a no-hair theorem appropriate to the supergravity theory involving SD-branes could 
be proven, however, the above solutions we have reviewed would be ruled out as candidates 
because their singularities are uncontrollably nasty. So we have to look elsewhere. 
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Let us take a brief sidetrip here to comment on the singularity story for supergravity 
fields of A^^l regular D-branes with timelike worldvolume. Certainly, the geometry of BPS 
D3-branes is nonsingular, and there are several other pretty situations known in the literature 
where branes "melt" into fluxes: the sources are no longer needed. However, the disappear- 
ance of D-brane sources for supergravity fields only occurs when the branes are BPS. If any 
energy density above BPS is added to these systems, singularities reappear: this certainly 
happens for the D3-brane system. Also, in the low-energy approximation to string theory, it 
is misleading to think of supergravity fields of Dp-branes as simple condensates of massless 
closed string modes. The reason is directly analogous to the fact that the Coulomb field 
of an electron cannot be a photon condensate because the photon is transverse. Similarly, 
Coulomb Ramond-Ramond fields of Dp-branes cannot be represented by supergravity fields 
alone. ^ This river runs deeper: in the decoupling limit, resolution of dilaton and curvature 
singularities for p 7^ 3 Dp-branes is in fact provided by the gauge theory on the D-branes 



24]. 



Let us now get back to our SD-branes. The supergravity situation looks similar to that 
for non-BPS (ordinary) D-branes: it seems that brane modes will be required for singularity 
resolution. Therefore, we are motivated to try to solve all problems with prior candidate 
SD-brane spacetimes by solving the coupled system of brane tachyon plus bulk supergravity 
fields with full backreaction. 

2.2 Unstable brane probes in sourceless SD-brane backgrounds 

The first progress towards the goal of singularity resolution in SD-brane systems was made 
by Buchel, Langfelder and Walcher [Q]. We now briefly review what is, for our purposes, the 
most relevant point of their work. 

Essentially, they take the reasonable point of view that the process of creation and sub- 
sequent decay of a SDp-brane must be driven by a single open string mode: the tachyon 
field denoted T{t), which lives on the associated unstable D(p+l)-brane. They use the p = 



non-dilatonic version of the worldvolume action |25, |27|, |l 



S brane 

J dP+^ye-^V{T)^-det {VG^p + d^TdpT) + /ip+i j f{T)dT A C^p+^\ 

(2.6) 

to study the dynamics of the tachyon and its couplings with bulk (closed string) modes. The 
operation V is for pullback, and a, (3 = 1 . . . (p+l). In section we will comment both on 
the validity of this type of effective action, and on the expected form for the tachyon potential 
V{T) and the Ramond-Ramond coupling f{T). For now, we just use it. 

The way we look at the calculation of ref. is as follows. Supergravity SD-brane fields 
should be regarded as arising directly from a large number of unstable branes. Then, using the 



^This is the case even though Dp-branes are "sohtonic" in string theory while electrons are fundamental in 
QED. The straightforward argument we use here depends only on the couplings of the charge-carrying objects 
to the bulk gauge fields. We thank Abhay Ashtekar for a discussion on this issue. 
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intuition gained from studying the enhangon mechanism [|30[ , it is natural to use an unstable 
brane probe to study more substantively the candidate supergravity solutions of refs. ^. 
Then, we look for problems arising in the probe calculation. The idea is that whenever the 
probe analysis goes wrong, it signals a pathology for the gravitational background. There 
are at least two ways the probe analysis can signal a problem: infinite energy or pressure 
density for the tachyon may be induced, {pprobe,Pprobe ioo), or there might not exist any 
reasonable solutions for T(t) across the horizon. 

So let us consider inserting an unstable brane probe in a background with fields corre- 
sponding to the sourceless SD-brane supergravity fields, eqs. {2A), of the previous subsection. 
The equation of motion for the open string tachyon is, generally, 

/ xi/ \£/ \i -<s>dV(T) ^ 1 

{-9tt)H9y^y^)n9y^y^)^e A2 



^ ^ d f ^^^^^ {9y^yl)H9y^y^)j \ ^ ^^.7) 



dt 



2 



where ' = d/dt, our notation for the Ramond-Ramond field strength is f(p+2) = dC^'P^^\ we 
have factored out Tp+i by including a factor of gs in C^p~^^\ and we also defined the following 
expression, 

A = l + ^. (2.8) 

9tt 

The question needing attention here is whether or not the tachyon field, regarded as a 
probe,^ is well-behaved when inserted in the candidate supergravity backgounds (|2.l|) . Ref. [@] 
provided a clear answer for the case of a d = 4 SDO-brane with dilaton field set to zero.^ Let 
us now see how this goes. 

The d = 4 SDO-brane background introduced in ref. pi] has the form. 



9tt — 9 — 0,0 2 ' 9y'y^ — U, 9xx — o C , l,^- 



and 

F2 = Qe2, $ = 0. (2.10) 

This spacetime metric has a regular horizon (a coordinate singularity) at t = uj, and a genuine 
timelike curvature singularity at t = 0. The expressions for the energy density (T/) and the 
pressure {T\ ) associated with the probe are respectively. 



Pprobe ~ , Pprobe ~ V{T)A^/^. (2.11) 



''The unstable brane will be a probe as long as its backreaction is small and can therefore be treated 
self-consistently as a perturbation. 

^The authors of ref. |Q| also investigated the effect of tachyon backreaction, but the ansatz they used for the 
supergravity fields was not general enough to handle our cases of interest. In particular, our general equations 
do not reduce to theirs upon consistent truncation. Also, their exposition of their backreaction analysis was 
extremely brief. 
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The only time (apart from t = 0) when the probe hmit becomes ill-defined is around 
t = uj. In fact, in the near horizon limit, the dynamical quantity A satisfies the simple 
ordinary differential equation 



This has the general solution 



where g is a constant of integration. There are two possible solutions at t = to: A = (g ^ 0) 
or A = 1 (g = 0). Clearly, the case for which A = corresponds to the probe limit breaking 
down since the energy density of the unstable brane diverges. The other possibility, A = 1, 
implies that the time-derivative of the tachyon diverges on the horizon. This last case is 
clearly pathological and cannot correspond to a physically relevant tachyon field solution. 

The conclusion is that unstable brane probes are not well-defined in the SDO-brane back- 
ground. Not only are they useless to resolve the timelike singularity at t = but, worse, 
they appear to generically induce a spacelike curvature singularity on the horizon at t = uj. 
That is, if we take the probe story to be a good indicator of the story for the full backreacted 
problem. 

One of our first motivations for the work leading to our paper was to plumb how restricted 
the conclusions of ref. |^ were. Did this above story work only for bulk couplings of the kind 
arising for SDO-branes in d = 4, where no dilaton field appears? Was it true only for the 
case of SDO-branes, which are a special case for SDp-branes since there can be no anistropy 
in a one-dimensional worldvolume? Are all timelike clothed singularities turned into naked 
spacelike singularities by probes? Could we even trust the probe approximation to tell us 
anything about the solution with full backreaction? 

The first generalization we considered was to look at an unstable brane probe in the 
background of the isotropic SDp-brane solutions of ref. However, what we saw there 
was that the naked spacelike singularities remained naked spacelike singularities; the small 
effect of the probe could not undo that pathology. Next, we moved to analyzing anisotropic 
solutions of the form those with regular horizons. Some of these are actually completely 
nonsingular; we analyzed the details of the probe computation in those backgrounds, and the 
specifics are recorded in Appendix^. The results there are simple to summarize: the solutions 
with singularities hidden behind horizons do not give rise to conclusions qualitatively different 
than what we have reviewed here for d = 4 SDO-branes. The picture therefore remains 
unsatisfactory. 

The upshot, then, is that the probe story does not resolve singularities found for sourceless 
supergravity SD-brane solutions. So we now move to the full backreacted problem for SDp- 
branes in d=10, which is the main content of our paper. 

3. Supergravity SD-branes with a tachyon source 

In this section we find, in the context of supergravity, the equations of motion associated with 




(2.12) 
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the real-time (formation and) decay of a clump of unstable D-branes. We begin this section 
by writing the form of the action which we will use in our analysis. We will concentrate on 
only the most relevant^ modes in both the open and closed string sectors; in other words, we 
keep in our analysis only the (homogeneous) tachyon and massless bulk supergravity fields. 
Potentially, the effect of massive modes could be encoded in a modified equation of state 
for the tachyon fluid on the unstable D-brane.'' We leave for future work the issue of non- 
homogeneous tachyonic modes, and of massive string modes in both the open and closed 
string sectors, for the coupled bulk-brane system with full backreaction. In order to use the 
supergravity approximation here self-consistently, we will take gg small but QsN large, and 
time-derivatives will be small compared to Ig ■ We will see that it is simple to choose boundary 
conditions in our numerical integration such that these remain true for all time. 

3.1 Preliminaries: action and equations of motion 

For this section we will be able to suppress R-R Chern-Simons terms in writing the bulk action. 
This is a consistent truncation, to set the NS-NS two-form to zero throughout the evolution 
of the system of interest, as long as consistency conditions on the R-R fields are satisfied. E.g. 
for the SD2-brane system with R-R field C^^^ activated, it is necessary to make certain that 
dC^^^ A dC^^^ = in order not to activate the NS-NS two-form and the accompanying Chern- 
Simons terms. Other cases are related to this one by T-duality. Therefore, we allow only 
electric-type coupling of the SDp-brane to C^^"^^^ (or equivalently magnetic-type coupling 
to CC^-P)). Later we will show that this ansatz is physically consistent provided we assume 
that there is ISO{p + 1) symmetry along the worldvolume of the SDp-brane, the object we 
are interested in. This is equivalent to considering only the lowest-mass tachyon, i.e., not 
allowing any excitations of the brane tachyon along the spatial worldvolume directions. Of 
course, the R-R field strengths are then very simple: G^^"*"^-* = dC^^^\ and the string frame 
bulk action takes the form [Kl^l 



L_ I d'W^gl^e-'" [n + 4{d^f] - \ } , (3.1) 



where IZ is the Ricci scalar. We use a mostly plus signature. In the above conventions, the 
R-R field solutions automatically get a factor of gg (as we mentioned also in the previous 
subsection), and 

le^Gio = (2vr)7<7?^8 , ^Dr> = —-\-^,. (3.2) 

Our analytical and numerical results in following sections will be given in string frame. If 
desired, it is easy to convert to Einstein frame - with metric g^u ~ with canonical normalization 
of the metric and positive dilaton kinetic energy, by using the standard d = 10 conformal 
transformation 

9_,. = e-^^-^-^/'g,.. (3.3) 

^Relevant in the technical sense. 

'^We thank Andy Strominger for this suggestion. 
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Stress-tensors are defined in Einstein frame, 



T = 



-1 5S, 



matter 



-g 



with the usual 



,(3-p)*/2 



G 



A2...A 



■p+2 



G 



2(p + 2) 



We can transform to string-frame "Einstein" equation using standard formulae 



(3.4) 

(3.5) 
(3.6) 

(3.7) 



For all matter fields except the dilaton, it is therefore obvious that string frame "stress- 
tensors" take the form 

rri 1 '^'S'matter g-^ 



'-HI/ 



^ 5g 



HI/ 



For the dilaton we find 



(3.9) 



and, obviously, familiar energy conditions for bulk fields are only satisfied in Einstein frame, 
not string frame. 

For the bulk field equations, we must include coupling to the brane tachyon — this is 
of course an important point of our paper. Hence, we now turn to the brane action. The 
brane theory is that appropriate to unstable D(p + l)-branes, for SDp-branes. We consider 
N branes. At low energy, the overall U{1) center-of-mass tachyon T couples as follows: 



A 



'brane 



|y dtdP+^y [- e-*V^F(r)] + J f{T)dTAGP+^^ , (3.10) 



IGttGio 

where the matrix Aa/s is defined as 



(3.11) 



where V stands for puUback. For the constants, our conventions are that T is normalized like 
Fafs, and also 

A = ^^^(IGttGio) = (iV5.)(27r4)'-^. (3.12) 

gs 

Notice that A is proportional to ggN. This will be the sole continuous^ control parameter 
associated with the physics of our final solutions for the coupled tachyon-supergravity system. 

*A is effectively continuous in the supergravity approximation since QsN is large and all derivatives small 
in string (is) units. 
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It is important to know when we can expect to trust the action we use. Our approach 
consists in assuming that the kinetic terms of the open string tachyon field are re-summed to 
take a Born-Infeld-hke form. Strictly speaking, this has only be shown to be a valid claim 



late in the tachyon evolution. We refer the reader to ref. [29| for more details on the limits in 
which this approximation holds (see also ref. ||3^). The functions V{T) and f{T) are therefore 
not known exactly at all times. For definiteness in our numerical analysis, we will choose a 



specific form and assume that the dynamics of the tachyon is governed by eq. (3.10). For a 



SD-brane, we make the choice V{T) = 1/ cosh(T/-v/2) = |/(^)| which has been shown to be 
the correct large-T behavior of the couplings. Our results turn out to be quite robust, in that 
their qualitative features do not depend on the precise form of V(T) and f(T). 

For the remainder of our discussion it will be convenient to use static gauge, which is an 
appropriate gauge choice for our problem of interest. Therefore, in the following, we will be 
rather cavalier about dropping pullback signs. 

When we get to solving the coupled brane-bulk equations, it will be convenient to allow 
for a density of branes, denoted p±, in the transverse space: 

/ dtd^+^y^ I dtdP+^xd^-Pxp±. (3.13) 

J brane J bulk 

Therefore our brane action becomes 

5brane = J d^'xp^ |_y(r)^e-* + f {T)e^'-^^^' q^^,„xp+2dx,]T} , (3.14) 

where e is the worldvolume permutation tensor with values (0, ±1). 

We are now ready to write down the coupled field equations. The simplest bulk equation 
to pick off is the dilaton. In string frame we see immediately that 

72^ /o^n2 , _ 1 Mp^A 



- {d<^r + 4^ = 8 (^^= ) yiT)V-A, (3.15) 



and for the Ramond-Ramond field 



(3.16) 



For the metric equation of motion (string frame "Einstein" equations), it is convenient to 
define 

7^'^, = n = n - -T^ (3.17) 

Therefore, we have 

Tim = -2^^^,$ - ^5^V2ci> + ^{d'^fg'^,. (3.18) 

For the brane stress-tensor, we need to figure out the dependence of V--A on g^i^ (the Wess- 
Zumino term clearly does not contribute). We find 

= -^e* (^) V{T)V^{A-'r^g^, , (3.19) 



2 vv^y 
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1 



-9 J 



1 



o 

The other object we need for the metric equation of motion is 



In Act, 



n[c] 



i/X...a 



(p + 2)! 



9Xa 



(3.20) 



(3.21) 



2(p + l)! ■^'^-^ 16 

Lastly, for the tachyon we find the equation of motion 

ViT)e-'^V^{A-^)^"'daT] + /(r)e^-^G^...A = . (3.22) 
In the Discussion section we will make some remarks about the robustness of these equations. 
3.2 The homogeneous brane self-consistent ansatz 

As pointed out earlier, we are interested in time-dependent processes by which massless 
Type Ila or Type lib supergravity fields are sourced by an open-string tachyon mode on 
the worldvolume of an unstable brane. A reasonable assumption is that the gravitational 
background generated by backreaction of the rolling tachyon is of the form 



(3.23) 



where the n-dimensional Euclidean metric dTj'^{k) is 

' dal forA: = +l 



dTi{k) 



dEl for A: = 
dHl for k = -l , 



(3.24) 



where is the unit metric on S*", the flat Euclidean metric, and dH"^ the 'unit metric' 
on n-dimensional hyperbolic space . The corresponding symmetry groups are 



SO{n + 1) for k = +l 
ISO{n) for yfc = 
SO{l,n) for/c = -l 



(3.25) 



1 and k\\ = 0, appears, after using an 



For /c= lb 1 we obviously require that n > 2. 

We note that the ansatz ( |3.23| ), for k_\_ = 
appropriate change of coordinates, to be equivalent to that considered for supergravity SD- 
branes in ref. |5|. One can show that this is actually not the case. In order to bring solutions 
of the form ( |3.23| ) with — oo < t < +oo to the form introduced in ref. |^] we must find a 
change of coordinates such that 



dt' 



dr^ F{T)'/\p{T)a{T))^p 

\a{r) 



(fci+fc)(p-l)/(7-p) 



(3.26) 



where — oo < r < +oo. It turns out that for all values of the parameters associated with 
the isotropic supergravity solutions, such a change of coordinates does not exist. However, 
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the main feature of our analysis is that we are ahowing for modifications of SD-brane physics 
in the region close to the spacelike worldvolume (the region around t = for the system of 
coordinates we use). It should therefore have been expected that our new solutions are not 
included in those presented in ref. ||5|. However, we do expect the asymptotics to agree. 

To be physically relevant, solutions should be asymptotically flat. For example, SD-brane 
gravity solutions will be such that 

lim a{t) = , lim R{t) = 1 , (3.27) 

t— >itoo t— »itoo 

for = and k± = —1. We will also see in our numerical analysis that only some values of 
k± and /cy are allowed. Also, we expect that for the dilaton and R-R field 

lim C{t)=0, lim 6(t) = 0. (3.28) 

+±oo t— >itoo 



By inspection of the tachyon equation of motion ( 3.22| ), we see that the electric- (or magnetic- 



) only ansatz referred to at the beginning of this section will be obviously consistent if we 
only allow worldvolume time-derivatives. This is tantamount to imposing an ISO{p + 1) 
symmetry on the worldvolume. Ref. (3^] argues that spatial inhomogeneities of the tachyon 
field will play an important role in the decay (a view which is also supported, although using 
a different line of reasoning, by the results of ref. |l^). It will be interesting to investigate 
the full importance of such effects in the context of our effective supergravity analysis. We 
will include a discussion of the nontrivial issues raised in the Discussion section. 

It turns out that the equations for the combined bulk-brane evolution in the time- 
dependent system are complicated enough to require numerical solution. For this reason, 
we will not be able to accommodate the most natural ansatz^ p± = 6{x). Instead, we will 
use the "smeared" ansatz also used by Buchel et al. in ref. Q, 

P± = Po^/g±- (3.29) 

Note that in this ansatz, the implicit time dependence in the transverse metric components 
is not varied in producing the equation of motion, rather it is only taken into consideration 
at the end of the calculation. Also, a smeared brane source does not contribute stress-energy 
perpendicular to the worldvolume, which is in the directions t ,y. It should be noted that the 
effect of using this ansatz will be minimized by using a small value for the density parameter 
Pq. Of course, the aim when using such an ansatz is to get rid of any brane action dependence 
on the transverse coordinates x. 

We should remark that supergravity solutions corresponding to unstable D-brane systems 
have been found before [^5|. Their solutions are time- independent, a feature which might seem 
rather unreasonable since they are, after all, supposed to describe unstable objects. Typically, 
these solutions are nakedly singular; there is no horizon. For reasons discussed previously. 



^Strictly speaking, instead of being a delta-function distribution, the more general ansatz for the source 
should be extended {e.g. a Gaussian) with its size of the order of the string length. 
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these solutions would therefore justifiably be regarded with some level of suspicion. Possibly, 
we should really regard these solutions as fixed-time snap-shots of the unstable brane system 
during its evolution. They do however reflect one desirable feature: taking into account 
warping of space in the directions transverse to the unstable branes. 

What we really want is a sort of hybrid of that approach - where transverse dependence 
is the only dependence ~ and what we are doing here - where time dependence is all there is. 
This is something we postpone to a future investigation; remarks on this will be given in the 
Discussion section. 

Let us now get back to the simplified ansatz, and just go ahead and solve it. We are 
therefore interested in the precise system of ordinary differential equations for our coupled 
tachyon-supergravity system. Using the form Ci2...p+i = C{t) (which is consistent with our 
ansatz) the equation of motion for the R-R field ( 3.16| ) becomes 



C + C 



XaF+^f{T)T. 



(3.30) 



Now let us find the dilaton equation of motion. A useful identity is 

7^ = 5(a$)2 - + ^^e'^G^ + \ (^) e*y(r) V^(A-i)N,. 

with which the dilaton equation of motion can be written, 
{P ~ 3) 2*^2 



4(p + 2)! ^ ' 



gxa 



(3.31) 



(3.32) 



This last expression is simply the Einstein frame equation of motion. So the dilaton in our 
ansatz satisfies 



$ + $ 



. R , ^ , a 
It a 



2$^ 



(3-p) / e^^C 



;.(P+l) 



(3-p)VA 



(3.33) 



This will be used whenever double time-derivatives of the dilaton need to be substituted for. 
With T = T{t) we find that the dynamics of the tachyon field is governed by 



T = Al^T -T 



{P+^)- 

a 



1 '^n^)+Zg)^^*^-(,+i)^^^ (334) 



V{T) dT V{T) 



where A = 1 — T^. In this paper we will be assuming that |/(T)| = V{T), a statement which 
has been shown to be correct only at past and future asymptopia. However, we have also 
done numerical experiments which show that some breaking of this relation at intermediate 
times (near the hilltop) does not change the important features of our solutions. 
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We now turn to the equations of motion for the metric components a{t) and R{t). For 
the stress-tensors, ehminating second order derivatives in matter fields, we have 



= 4($)2 - 2$ 



p+1 - + (8-p)- 
a H 



+ 



(5 - 2p) / e^C 



+\xe''V{T) 



{7-2p)VA- 



y a 4 I aP+i 



lxe''V{T)VA, 



R 4\aP^ 



A +lAe*y(T)^/A. 



(3.35) 



The components of the Ricci tensor are easily evaluated: 



7^* 



7^^ 



/ X o / ^ R 
(p+l)- + (8-p)-, 
a It 



d . .aR 
- + (8-p)-- +p 
a a R 



+ 



R , _.dR ,„ , 
R^^P^'^aR^^'-P^ 



R\ _fc± 
R R^ 



(3.36) 
(3.37) 

(3.38) 



For the (string-frame) "Einstein" equations, the time, longitudinal and transverse components 
are respectively 



(P + 1)- + 
a 



R 

P^R 



(P+1)- + 



R 

P^R 



i7-2p)VA 



,(3.39) 



a 
a 



-{8-p) 



aR 
aR 



p 



2 ku 




^Ae*y(T)^/A, (3.40) 



R 
R 



-iP + l) 



Ra 
Ra 



i7-p) 



+ 



i?2 



+ 2^>— -h - 
R 4 



+ ^Ae*y(r)VA. (3.41) 



In the end we have a system of five second order coupled ordinary differential equations 
for {T,C,^,a,R} as functions of t. These are respectively equations ( [3.34 ), ( 3.3C| ), ( ^.33 ), 
( |3.4[l| ) and ( 3.41]) . For consistency this system of equations must be supplemented with the 



- 15 - 



first-order constraint 



2(p+l)(8-p)-^+p(p+l) 
a R 



a J 



-40 



(p + l)- + (8-p)- 



+ {7-p){8-p) 



R 



+ 



i?2 



(3.42) 

obtained by plugging eqs. ( p. 40] ) and (3.41) in eq. (3.39). Of course, if this last equation is 
satisfied at t = it will be for all times. 



4. The roll of the tachyon 

We refer to a supergravity SD-brane as the field configuration generated by a system composed 



of a large number, N, of microscopic SD-branes. As emphasized in section 3.2, there is 
a single continuous parameter that controls the dynamics of these fields, i.e., A = po^- 
This is the parameter that determines the relative importance of the unstable brane source. 
Clearly, for A — > the open string tachyon decouples from the bulk fields (no backreaction) 
and the corresponding supergravity solutions will presumably be the singular ones found in 
refs. 

In this section we present solutions with the symmetries of SD-branes and a non- vanishing 
A. We demonstrate that they are generically non-singular. 

4.1 Tachyon evolution in flat space 

The solutions relevant for SD-brane physics in Type II a,b superstring theory could be the 
ones corresponding to an open string tachyon evolving from 

lim T{t) = +CX) (4.1) 

t— >— oo 



to 



lim T{t) = -oo (4.2) 



in a symmetric runaway potential of the form shown on figure |T[ Solutions of this type 
correspond to the tachyon evolving between two different minima of the potential. Possible 
initial conditions (at t = 0) for these solutions are of the form 

r(0) = const. , r(0) = . (4.3) 

Another set of solutions corresponds to a tachyon evolution with 

lim T{t) = aoo, (4.4) 

t— utoo 
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where a = 1 and a = — 1 lead to equivalent solutions. Appropriate boundary conditions for 
the tachyon (at t = 0) are then of the form 



r(0) = 0, r(0) = const. 



(4.5) 



In this section we characterize the supergravity solutions generated by a tachyonic source 
with the properties mentioned above. The results we present are for solutions with boundary 
conditions (4^) and asymptotic behavior ([4.4|). We have found that the solutions associated 
with the boundary conditions (|4.3|) and asymptotic behavior ([4.1|) have similar qualitative 



features. Strictly speaking, our approach to studying real-time evolution in supergravity can 
also be extended to more general cases, i.e., brane decay or creation (half SD-branes). 



/ 0.6- 




y(T) : 




/ 0.4- 




/ 0.2- 





-A 



Figure 1: Possible form of the open string tachyon potential V{T) on a SD-brane in Type II a,b 
superstring theory. 



In our analysis we use the potential 

V{T) 



1 



(4.6) 



cosh (r/^/2) 

because it agrees with open string field theory calculations for large values of the tachyon for 
unstable D-brane systems in Type II a,b superstring theory.^'' It is not known what the exact 
potential is for intermediate times but our solutions are only mildly affected by its particular 
form. The expression for the tachyon when A = 0, i.e., when there are no couplings to the 
bulk supergravity modes, is^^ 



T(t) = —V2 arcsinh ( — T(0) sinh 



V2 



(4.7) 



^°In bosonic string theory the potential is asymmetric and unbounded from below as T — > —00. 
^^We refer the reader to Appendix H for a derivation of this expression. 



-17- 



for the boundary conditions ([4.3| ). When the boundary conditions are of the form ([4.5| ), the 
analytic expression for the tachyon is 



T{t) = — \/2 arcsinh | sinh 



■ — I cosh —pz 
V2 J V2 



(4. 



Figure |2| shows a tachyon profile for T(0) = 0.83 and T(0) = 0. Homogeneous solutions 
such as eqs. (4.7) and ( |4.8| ), derived from a tachyonic action, correspond to a fluid that has 
a constant energy density and vanishing pressure asymptotically (tachyon matter). We will 
shortly see how these features are affected by the inclusion of couplings to bulk modes. 



\ 




\ ^ 




\ T(t)6- 













-10 



10 



Figure 2: The tachyon field evolution for r(0) and r(0) = 0.83. 



4.2 R-symmetry group 

Before presenting the numerical results we comment on the issue of R-symmetry. As pointed 
out in ref. SD-brane solutions should be invariant under the transverse Lorentz trans- 
formations leaving the location of the brane fixed. This corresponds to an 50(1,8 — p) 
R-symmetry. This property is embodied in the supergravity solutions found in refs. [||, ^, ^ 
where the transverse space metric has a factor of the form: R(t)'^dHg_p. The embedding 
group of the hyperbolic space H^'^ is S0{1, 8 — This explains why we study supergrav- 
ity solutions with fcy = and k± = — 1 in more details. However, the cases with k± = are 

^■^The interpretation in terms of R-symmetry is relevant for the idea that SD-brane gravity fields are holo- 
graphically dual to a worldvolume Euclidean field theory. 

^^The intuition for the nature of the R-symmetry group is inherited from the AdS/CFT correspondence. 
For example, the metric of a 3-brane is of the form: ds^ = /(r) [-dt^ + d-if] + l//(r) [dr^ -I- r^dQ.^] . The 
R-symmetry group in this case is SO{6), a statement that can be traced to the fact that the near horizon limit 
of the geometry is AdSs x S^, i.e., the gravitational background dual to the worldvolume field theory on an 
ensemble of D3-branes. 



- 18 - 



also candidate solutions for SD-brane (and, more generally, unstable D-brane) supergravity 
solutions. We present an analysis of those and other cases in section [4.4| . 

It is suggested in ref. [P that the spacelike naked singularities |jl], ^, |5) associated with the 
supergravity SD-branes could be resolved by using a metric ansatz that allows for a breaking 
of the R-symmetry in the region around the core of the object (t = 0). The intuition from 



the AdS/CFT correspondence comes from ref. |36| which describes cases of spontaneously 
broken R-symmetry. Our ansatz as it is cannot accommodate such an R-symmetry breaking. 
Presumably, this would correspond to a time-dependent process by which the R-symmetry 
is broken down to SO {8 — p) in a finite region. The closest our ansatz can come to realizing 
this scenario is if we consider the cases with k± = 0. Then, the transverse symmetry group is 
ISO{8 — p), the compactification of which is S0{8 — p). We find that solutions with k± = 
are regular and asymptotically flat. Moreover, as we will see in section [4.4| , a generic feature of 
the k± = solutions is that the effective string coupling exp($(t)) vanishes asymptotically. 
We will also remark on the cases with k± and/or fcy equal to +1, which all develop big-crunch 
singularities in finite time. 

4.3 Numerical SD-brane solutions 

Our ansatz for describing supergravity SDp-brane solutions is 

dsgp = -dt^ + a{tfdf + R{tfdHi_p, 

C = Cit), <D = cI>(t), (4.9) 

where we have taken fcy = and A;_l = — 1 in accord with the original paper ref. Again, 
these supergravity bulk modes will be excited by couplings to an homogeneous open string 



tachyon as described in section |3.1| . We consider only homogeneous solutions (i.e., depending 
only on time). As discussed earlier, it is possible that inhomogeneities might play an important 
role in the creation/ decay process |3^, but we are postponing this issue for now. 

The resulting system of coupled differential equations for which we will find solutions is 
of course highly non-linear. Among other things, this implies that it is not possible to extract 
scaling behavior for the fields from their equations of motion. We could therefore expect that 
the behavior of these solutions depends in a physically important way on the initial conditions 
for the field components involved: a{t), R{t), C{t), <I>(t) and the source T{t). We will see 
that the qualitative behavior of the supergravity SD-brane solutions actually does not depend 
very much on these initial conditions. 

The multiplicity of potential solutions could be large: one solution would be expected for 
each set of initial conditions. Now, for SD-branes the number of solutions is greatly restricted 



by the constraint equation (|3.42D . It would appear natural to consider, as candidate solutions 
for SD-branes, those which are time-reversal symmetric around t = O}^ In fact, precisely 



14 



We thank Alex Buchel for pointing out that the time-reversal symmetric solutions are inconsistent. In 



fact, the constraint equation (3.42) cannot be satisfied (for the case fcy = and k± = —I) with the boundary 



conditions: d(0) = 0, i?(0) = 0, (p{0) = 0, C(0) = and r(0) = 0. 
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time-symmetric solutions actually do not exist! That is, unless they have /cy = 0,k± = +1, in 
which case they have R-symmetry which is markedly in conflict with the proposals of earlier 
papers, e.g. As we said above, these solutions with wrong R-symmetry also develop a 
big-crunch singularity in finite time. 

In a sense, motivated by inflationary cosmology, the absence of precisely time-symmetric 
solutions is not particularly bothersome to us. What we mean by this is that it is an extremely 
fine-tuned situation to have exactly zero kinetic energy in each of the bulk fields at t = 0. Any 
small quantum fluctuation of a bulk field takes us away from time-symmetry. Therefore, all 
solutions which we will exhibit here will have some kinetic energy in one or more of the bulk 
fields at t = 0. And the kick in bulk field(s) at i = 0, required to solve the initial constraint, 
can actually be made very small by, e.g., choosing R{0) large. Such initial conditions are 
quite generic: even a little bit of kick in a{t) alone will suffice to give completely nonsingular 
solutions for the entire time-evolution. The precisely time-symmetric solutions are however 
inaccessible; see subsection 4.4.2 for a more detailed exposition of this case. 

Another reason why we find a slight bulk asymmetry at t = reasonable relates to 
particle/string production. We have of course neglected such production in our analysis. It is 
nonetheless clear that, for a full SD-brane solution, backreaction will combine with particle 
production to make bulk fields naturally asymmetric. It is only in the approximation of zero 
backreaction that time-symmetry is possible when particle production occurs, but of course 
that approximation cannot be self-consistent. 

Let us now move to the solution of our problem of interest. Unfortunately we have been 
unsuccessful at finding analytical expressions for the bulk modes and tachyon field when their 
mutual coupling is non- vanishing (A 7^ 0). We therefore resorted to solving the corresponding 
system of differential equations numerically.^^ In what follows we show the results associated 
with a SD4-brane, and present some interesting analysis of the effect of varying the initial 
conditions. We provide general comments for other SD-branes with p < 7 and explain how 
the SDp-branes with p = 7,8 are different. We comment on the robustness of the solutions to 
changes in initial conditions. Throughout our analysis we pay special attention to the impact 
of varying the initial condition <&(0) on the dilaton, because this controls the string coupling 
close to the hilltop. 

For finding the numerical solutions we fix the initial conditions at t = 0, and evolve this 
data forward in time. Then, we evolve the same data backward in time from t = 0. The 
result is the solution associated with a full SD-brane, i.e., the bulk fields sourced by the open 
string tachyon as it evolves from t = — co to t = +00. In this section we present the results 
associated with a SD4-brane with boundary conditions 

r(0)=0.83, r(0)=0, a(0) = l, d(0)= 0.091, i?(0) = 10, 
i?(0)=0.01, (f){0) = -l, 0(0) = 0.01, C(0) = 0, (7(0) =0.01, (4.10) 

and A = 0.1 . Again, these solutions correspond to the tachyon rolling up the potential from 
T = +00 {t = —00) and coming to a halt for T(0) = 0.83 which corresponds to a turning 
^^Recall that this is the primary reason why it is so difficult to include ^-dependence in our ansatz. 
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point. Then, the tachyon evolves toward the bottom of the potential at T = +00 for t = +00. 
This type of tachyon evolution was considered, for example, in refs. [14, 10| where they are 
called full S-branes. We could of course have chosen different initial conditions - including 
some with less time-asymmetry; these are used just for the sake of illustration. 



4.3.1 Deformation of the tachyon field 

Figure || shows the evolution of the tachyon for the S4-brane. For A = (no coupling between 
the closed and the open string modes), we find 

lim T{t) = ±t. (4.11) 

t— >itoo 

This observation is suggestive that the tachyon field might play the role of time itself in 
cosmological models driven by brane decay (as proposed in ref. ||2^). We find that this 
asymptotic behavior for the tachyon survives when couplings to bulk modes are introduced. 
In fact, for A 7^ we find 

lim T{t) =±{t + K±T)- (4.12) 

The constants k±t depend non-trivially on the boundary conditions at t = and p. For 
example, as p is increased Kt decreases. Also, for large values of <I>(0) the tachyon deformation 
from the flat space case becomes larger. As mentioned in Appendix]^, for A = the state of 
the tachyon for large time {t — > ±00) is that of a perfect fluid with constant energy density 
and vanishing pressure. This is the so-called tachyon matter. We find that for A 7^ 0, both the 
energy density and the pressure (physical quantities measured in the Einstein frame) vanish. 
In other words, the tachyon matter is clearly only an illusion of the Qs ^ limit. 

We consider briefly the effect of varying the initial conditions T(0) and T(0). For half 
SD-branes {i.e., the future of the full SD-branes) we find that the time it takes the tachyon 
to reach the bottom of its potential increases for smaller values of r(0) and r(0). Not only 
that, for very small initial velocities the tachyon stays perched at the top of the potential for a 
certain period of time. In general, we observe that it takes less time for the tachyon to reach 
the bottom of its potential when we increase A. Now, we also observe that the difference 
between the curves associated with flat space (A = 0) and A 7^ tachyons decreases as T(0) 
and r(0) increase. Also, for large negative values of $(0), Kt becomes very small. This is 
simply a reflection of the fact that such cases correspond to a very small initial string coupling 
(see below for details). 

There is another interesting feature of the tachyon when coupled to the massless closed 
string modes. Firstly, the time it takes the tachyon to reach the bottom of its potential is not 
significantly altered even when considering large values of the 'coupling' A. Finally, we find 
that as the coupling A is increased, Kj., or the deformation away from the flat space tachyon, 
increases correspondingly. 
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4.3.2 Time-dependent string coupling 

The string coupling is given by the expression 



Qs = e<*°> , (4.13) 

where $o is the background dilaton field in the absence of sources, i.e., strings and D-branes. 
Typically the presence of stringy excitations will modify the coupling of the theory, 

5s-5se*(''), (4.14) 

where is a spacelike variable for D-branes and is timelike for SD-branes. 

For supergravity Dp-brane solutions the dilaton field is (see, for example, ref. [p7|]), 

/ ,7-p\ 3(3-P) 

9s{r) = 9se''^^^=gsU + '''l7; j , (4.15) 

where Cp = {2^/^T)^^^^[^{7 — p)]- The string coupling is seen to vary according to whether 
test closed strings propagate close or far from the horizon (r = 0) of these geometries. For all 
static supergravity solutions (including NS5-branes), the asymptotic string coupling is such 
that 

lim gs{r) = g^. (4.16) 
r— >+oo 

The effect of supergravity D-branes is therefore to modify the coupling locally. For p < 3, it 
is large close to the horizon and decreases to gf^ as r — > +cxd. For p > 3, the string coupling is 
small close to the horizon but increases to gs for large r. The case p = 3 is special because the 
dilaton field sourced by the 3-brane is constant throughout the spacetime. Typically, the size 
of the region where the dilaton is not constant depends on the parameter gsN . Large values 
of this parameter are associated with larger regions where dilaton perturbations associated 
with the brane are noticeable. 

The solutions associated with supergravity SD-branes induce dilaton perturbations cor- 
responding to a time-dependent string coupling, 

5.(t)=ff.e*W. (4.17) 

We find that the time dependence of the dilaton sourced by SD-branes is qualitatively different 
when compared to the radial dependence of the dilaton associated with regular D-branes. 

Figure ^ shows the time-evolution of the dilaton function 57 ^e*^*^ for a SD4-brane with 
boundary conditions ( 4.10| ). Typically, the function g~^e^{t) decreases from t = as t — > ±00. 



We find that smaller values of p correspond to larger asymptotic string couplings. 



^^This is an other example where features of SD-branes are not simply those inherited by analytic continu- 
ation of D-branes. In fact, a double analytic continuation (r — > ir and t it) of the supergravity Dp-brane 
solutions lead to objects with an imaginary R-R charge. This is unphysical. 
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The dilaton field generated by SD-branes has at least two interesting properties. First, 
all solutions are such that the dilaton stabilizes to a constant asymptotically, 

lim $(t) = $±^. (4.18) 

More generally, the relation between $+00 and $-00 depends on the initial conditions. This 
last statement applies to all other bulk fields. Secondly, the asymptotic value of the dilaton 
is always smaller than its initial value at t = 0, 

^±00 < $(0). (4.19) 

This implies that the late/early string coupling {t ±00) is always smaller than the coupling 
when the tachyon is at the top of its potential (t = 0), i.e., 

g,e*±- <5.e*W. (4.20) 

We find that the qualitative features shown on figure ^ are preserved when the boundary 
conditions on the various fields are changed. Nevertheless, we consider the eff'ect of varying 
<&(0) in some detail. An interesting quantity to study is the ratio 

^ = (4-21) 

which gives a quantitative measure of how much the initial string coupling is modified asymp- 
totically. The tachyon profile is not altered significantly when the initial condition on <I>(t) 
is varied. Nevertheless, we observe that that for large values of $(0) (large initial coupling) 
the tachyon field reaches the bottom of the potential well faster. A large initial coupling also 
means that the bulk fields relax faster to their stable asymptotic configuration compared to 
cases where <1>(0) is smaller. The overall eff'ect on the bulk fields is also enhanced for larger 
values of $(0). For example, as the initial coupling is increased the scale factor o(t) stabilizes 
to significantly smaller values (see next section). As for the dilaton field itself, we find that 
the ratio h is large for larger values of $(0). For very small values of the initial string coupling 
(large negative values of <&(0)), the ratio h approaches unity. 

In summary, we find that the parameter <^(0), i.e., the parameter determining the string 
coupling when the tachyon is close to the top of its potential, strongly determines the impor- 
tance of the unstable brane source effect on the supergravity bulk modes. 

4.3.3 Gravitational field 

We now describe the effect of the unstable brane source on the time-dependent metric com- 
ponents a{t) and R{t). 

Figure |^ shows the scale factor a{t) on the worldvolume of a SD4-brane with the boundary 
conditions (4.10). A general feature is that away from t = the scale factor monotonically 



decreases from its initial value, a(0), to a stable asymptotic value, 



lim a{t) = a±oo, (4.22) 

>itoo 



''Obviously ^(O) can be taken to have negative values. 
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Figure 3: Time dependence of the function e*'*' / gg for a SD4-brane with boundary conditions ( 4.10 ) 
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Figure 4: The scale factor on the worldvolume of a supergravity SD4-brane with boundary conditions 



with a+oo > O-oo- The time it takes for this scale factor to reach its asymptotic value 
corresponds roughly to the time it takes for the tachyon to reach the bottom of its potential. 
As pointed out above, for large initial values of the dilaton, <I>(0), the asymptotic values of 
the scale factor, a±oo) become smaller. Correspondingly, when the initial coupling is weak 
the effect of the probe on the bulk modes is small and a{t) stabilizes to a value closer to its 
initial value a(0). 
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Figure |5| shows the behavior of the metric function R{t). A generic feature of the super- 
gravity SD-brane solutions is that 



hm R{t) = ±{t + 

4— >itoo 



(4.23) 



The constants k^^j are generically larger for larger values of p. 
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Figure 5: The SD4-brane transverse scale factor R{t) with boundary conditions (4.10) 



4.3.4 Ramond-Ramond field 

Figure ^ shows the time dependence of the Ramond-Ramond form field for a supergravity SD4- 
brane with boundary conditions ( 4.10| ). Again, the energy stored in this field (proportional 



to its time-derivative) goes to zero in approximately the time it takes for the tachyon to reach 
the bottom of its potential. A generic feature of the Ramond-Ramond field associated with 
a SD-brane is therefore, 

lim C(t) = C±oo, (4.24) 

t— >±oo 

where C±oo is a constant. Typically we find that these constants are smaller for larger values 
of p. 

4.3.5 Curvature bounds and asymptotic flatness 

The supergravity equations of motion are derived from a worldsheet calculation by requiring 
that, at a certain order in perturbation theory, the beta-functions associated with bulk fields 
vanish. Typically, there are higher order (in a') curvature corrections to the beta-functions. 
These corrections are negligible only if the curvature involved is small when measured with 
respect to the string length. Is = \fa' . This is why our solutions can, strictly speaking, be 
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Figure 6: The SD4-brane Ramond-Ramond field C{t) associated with the boundary conditions (4.10) 



trusted only if the curvature involved is such that \Tl\ , \TZ^uT^^^ \ , \T^t^upxT^^'^''^\ are small. 
We verify that this condition is satisfied by studying the behavior of the time-dependent Ricci 
scalar of the supergravity SD-branes, 

7^(^) = 2ip + 1)- + 2(8 - + 2(p + 1)(8 - 

a K a K 

where k± = —1 and fcy = for the cases of interest here. 

A property of the solutions, which is apparent from studying the evolution of bulk modes, 
is that of asymptotic flatness. In fact, we find 

lim dsl =-dt^+al^df + {t + K±RfdHl (4.26) 
hm c^(f) = $±oo, lim C(t) = C±oo, (4.27) 

£— >±00 t— >±00 

where a±oo, k±ii, $±oo and C±oo are constants. Both the first- and second-derivative of the 
bulk modes vanish asymptotically. The resulting brane configuration is then clearly flat for 
t — > iboo, as it should be. 

An important question to answer at this stage is: Which quantity in the problem sets 
an upper bound on the curvature for SD-branes ? First, we have found that whatever the 
curvature is at t = 0, its absolute value will never exceed it significantly in the course of 
the evolution. This is true only for half S-branes with boundary conditions corresponding to 
positive derivatives of the bulk fields. More generally for full S-branes the acceptable solutions 
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are those with a combination of the boundary conditions such that 



linl|7^(^)| small . (4.28) 

These requirements are certainly attainable within the self-consistent supergravity approxi- 
mation we are considering — along with the specific ansatz we introduced in order to be able 
to solve the equations numerically. 

4.3.6 The p = 7 and space-filling SD-branes 

The case p = 8 should be special because there is no transverse space into which "energy" 
can be dissipated. The only bulk fields involved are then the metric component, a{t), the 
dilaton, ^{t), and the Ramond-Ramond field, C{t). We find that the time-derivative of the 
scale factor, d(t), decreases to zero, as t — > ±00, many orders of magnitude slower than for 
p < 7. The scale factor a(t) also goes to zero after an infinite time (see section |4.3.7 for a 



physical interpretation) which is to be contrasted with the cases p < 7 where a±oo 7^ 0. One 
would expect that to be the source of a curvature singularity at t = ±00 but it is not the 
case. The Ricci scalar for the space filling SD-brane is 

nit) = -72 f - V - - f ^\ + 366- + -Ae* Ai/V(r). (4.29) 
\a J 2 \ I a 2 

Curvature singularities are avoided because all time-derivatives in the problem go to zero 
faster than the scale factor as t — > ±00. In particular, the quantity e^C goes to zero faster than 
a{t) for large time. Also the string coupling slowly goes to zero asymptotically ('I'ioo ~^ — cxd)! 
For the Ramond-Ramond field we find the same behavior as for the cases p < 7. The only 
difference is that the relaxation time of the bulk modes is many orders of magnitude larger 
than in the other cases. 

The functions a{t), ^{t) and C{t) associated with the SD7-brane behave in the same way 
as the space-filling SD-brane. Of course in this case there is a transverse space and we find 
that the relaxation of |-R(i)| to its asymptotic form t + k^r takes an infinite amount of time. 
In other words, it relaxes to its asymptotic form much slower than for the cases p < 7. 

4.3.7 Einstein frame 

Let us end this section with a remark about the physics of our SD-branes in Einstein frame. 
This would be the more natural and more physical frame to use in discussions of potential 
uses of rolling tachyons in the context of cosmology. 

The transformation from string frame to Einstein frame involves a multiplicative factor of 
e-*/2 for d = 10, which is the dimension in which we are working here. Now, we have already 



commented at length on the behavior of the time-dependent dilaton field in section 4.3.2. 
The essential physics there was that the dilaton is biggest at the top of the potential hill; in 
particular, it stabilizes in the infinite past and future at a constant value smaller than the 
initial condition at the hilltop. Dilaton derivatives also remain small at all times during the 
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evolution. Therefore, most of the quahtative features of our solutions will be preserved upon 
transformation to Einstein frame; in particular, all solutions remain completely nonsingular. 

A case that deserves further comments is that of the space filling SD8-brane discussed 
earlier. In string frame the metric is written 

dsls = -df + a{tfdf. (4.30) 
Upon converting to Einstein frame we get the metric 



where we have used a change of coordinate such that dr^ = e and where 



For the SD8-brane we found 



dsEss = -dr"^ + CLEi'rf'df', 



lim a{t) = , lim e*^*) = 0. 

t— >itoo t— >itoo 



(4.31) 
(4.32) 
(4.33) 



In string frame this implies that the metric "closes off" at infinity but in the (more physical) 
Einstein frame the converse happens, i.e., the limit r — > ±00 corresponds to a constant scale 
factor, 

lim a£;(r) = const. (4.34) 

This is a very sensible result. We mentioned before that, because of the absence of a transverse 
space, there appeared to be no channel into which the energy could go. We therefore find 
that asymptotically the energy has gone into inflating the worldvolume of the SD-brane. In 
fact, the relaxation time for the gravitational field is essentially infinite. 



4.4 Other classes of solutions 

In this section we present a more general analysis of the family of solutions with parameters 
{^11 ) k±} associated with the metric ansatz ( 3.23| ). 



4.4.1 Comment about the numerical analysis 

For A = there is clearly no source for the the bulk fields. The equation of motion for the 
tachyon is then ( [8.101 ). In Appendix ^ we found analytic expressions for the corresponding 
solutions. They have the large time property 

^^^^ = const. (4.35) 



However, if we solve eq. ( B.IOD numerically (using the same techniques as in section Q) we 



find that V{T)/\/A appears to become constant but, after some more time has elapsed, its 
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behavior becomes erratic, i.e., it starts oscillating with increasingly large amplitudes. Clearly 
this is only an artifact of the numerical analysis. The tachyon evolution is such that 



lim V{T) = 0, lim \/A = 0. (4.36) 

t— >±oo t^±oo 

The source of the problem can be traced down to the numerics having difficulties to evaluate 
a g division at large times. 

The same thing happens for A 7^ 0. The quantity that becomes ambiguous for large time 
is then 

(4.37) 

Again, the problem is associated with the fact that both e*^*-* '^{T) and -v/A are zero for large 
time. It is important to resolve this ambiguity because the quantity ( 4.37| ) directly feeds in 



the equations of motions for the bulk fields. We are able to show numerically that (at least 
for the solutions considered in this work) for large time 

^^ = 0. (4.38) 



Similarly to the A = case, past some finite time the behavior of the function ( |4.37 ) becomes 



erratic. To obtain numerical solutions representing the evolution for all times we did the 
following: past the time value where ( |4.37| ) becomes ill-behaved, we solve the system of 
differential equations without a source, i.e., for 

£!^1^^ = 0, e*(*V(r)VA = 0, (4.39) 

and appropriate boundary conditions. This introduces negligible errors. 
4.4.2 Time-reversal symmetric solutions 

We consider the time-reversal symmetric solutions, i.e., those associated with bulk fields 
having vanishing time-derivatives at t = 0: 

d(0) = i?(0) = C{0) = 6(0) = 0. (4.40) 

The constraint equation ( |3.42| ) at t = is then 

^ ^p(p+l)^ + (7-p)(8-p)|^V (4.41) 



A Ae* V^'^ 'a2 ^ ^^ii2 



The LHS being positive-definite, the constraint can be satisfied if and only if the metric 
ansatz contains a subspace of positive curvature. These consist in five categories of solutions, 
i.e., = {0, 1}, {— 1, 1}, {1, — 1}, {1, 1}, {1, 0}. We have performed a detailed analysis 

of these solutions for p = 4 and r(0) < 1/10. Typically, we find that when the tachyon 
has reached a point of its evolution where |T^(T)| ~ 0, labelled tc, the solutions develop a 



- 29 - 



curvature singularity. The behavior of the bulk fields is as follows. The time derivative of the 
scale factor is such that 

lim a = -oo, lim a{t) = 0, (4.42) 

t >tc t ^tc 

which corresponds to a big-crunch singularity on the {p + l)-dimensional worldvolume. The 
cases = —l,k± = 1}, {A;|| =0,k± = 1} and {k\\ = l,k± = 0} are such that 

lim R{t) = -oo, lim R{t) = 0, (4.43) 

corresponding to the transverse spherical space collapsing to zero-size in finite time. For 
{A;|| = l,k± = 1} and = l,k± = —1}, we find 

lim R{t) = +00, lim R{t) = +oo. (4.44) 

For k± = —1, R(t) goes to infinity faster than t. All solutions are such that 

lim 6(i) = -oo. (4.45) 

t >tc 

Generically, the large time behavior of the R-R field is well-behaved, i.e., 

lim C{t) « 0, lim C{t) = Co, (4.46) 

where Cc is finite but typically many orders of magnitudes larger than the constants Coo 
associated with the regular solutions presented earlier. The case = 0,k± = —1} is special 
because then the time derivative of the R-R field diverges as well at finite time. 

We believe our conclusions to be unaltered for other values of p. We found no evidence 
that the singularities are resolved when the time-reversal symmetry is broken. 

4.4.3 More regular solutions 

Another class of candidate SD-branc solutions we have studied are those with fey = and 
k± = 0. The results we present for these solutions are generic, i.e., they hold for all p and 
reasonable boundary conditions (i.e., first derivatives not too large). Firstly, the solutions are 
always asymmetric around t = since it was shown before that, for consistency, at least one 
of the bulk fields must have non- vanishing kinetic energy at t = 0. The level of asymmetry 
will be reduced by having smaller derivatives of the bulk fields at t = 0. The corresponding 
solutions are regular and asymptotically flat. We find that 



lim a(t) = 0, lim a{t) = 0, (4.47) 



and 



lim R{t) = 0, lim R{t) = const. (4.48) 

i— »±oo t— >±oo 

In Einstein frame both scale factors asymptote to non-zero constants. For the dilaton we 
obtain 

lim i{t) = 0, lim $(t) = -oo, (4.49) 

t— >±oo t— >±oo 
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i.e., the string coupling vanishes asymptoticaUy. FinaUy the R-R field behaves like 

lim C{t) = 0, lim C{t) = const 

t— >itoo t— »itoo 

We note that the relaxation time for these solutions is many orders of magnitude larger than 
for the {fcy = 0, k± = —1} cases presented earlier. 

The two remaining cases are {A;|| , k±} equal to { — 1, —1} and { — 1, 0}. We found evidence 
that the corresponding solutions are regular and asymptotically flat. 

5. Discussion 

Our primary motivation for this work was the general problem of seeking mechanisms for 
resolution of singularities in spacetimes of interest in string theory. SD-brane supergravity 
spacetimes presented in refs. [|l|, ^ ^ create somewhat of a supergravity emergency because 
they are not only singular but nakedly so. 

An important first step in the resolution program was made in ref . Q , where some effects 
of unstable brane probes in these backgrounds were considered. In particular, the probe 
physics was also sick, and taking this analysis seriously led to an even more dire assessment 
of the likelihood of singularity resolution without resorting to the inclusion of massive string 
modes in both the open and closed sectors. In our work, then, we began by plumbing the 
depths of the probe approach. We found it to be generally insufficient for our purposes; one 
reason is that the probe approximation takes itself out of its regime of self-consistency. 

We then launched into an investigation of the physics of the gravitational fields exerted by 
SDp-branes for general p by including backreaction. In order to get started on this problem, 
we had to make the approximation of considering only the most relevant open- and closed- 
string modes, with full gravitational backreaction taken into account. The equations we 
derived are highly nonlinear and couple brane with bulk, so did not lend themselves to solution 
analytically. We therefore resorted to numerical techniques to search for solutions. Because of 
this restriction, we had to use an ansatz which smeared the branes in the transverse space; this 
allowed us to turn the equations into ODE's and integrate them numerically. An essential step 
was to begin the numerical integration near the top of the potential hill, and then reconstruct 
asymptopia, which we were able to do successfully. Generically, the solutions are time-reversal 
asymmetric. We have shown that the time-reversal symmetric solutions with the correct R- 
symmetry are unavailable. Given our two-stage approximation (lowest modes, and smeared 
ansatz), we found it rather satisfying that significant progress in the resolution program is 
already found at this level. In particular, our solutions for rolling tachyons backreacting 
on spacetime are completely nonsingular, and our approximations satisfy the fundamental 
property of self-consistency. We find these conclusions suggestive of resolution of the SD- 
brane spacetime singularity emergency. 

It is however hard to know for sure whether our nonsingular results will survive re- 
finement. Therefore, let us now make some specific remarks about technical roadblocks we 
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encountered which forced us to make approximations, their physical consequences, and future 
outlook. 



Section 3.1 was where we derived the coupled tachyon-supergravity equations for a general 
brane distribution, assuming Chern-Simons terms are turned off in a consistent truncation. 
Our resulting equations are on the one hand remarkably non-robust, and on the other quite 
robust. What we mean by non-robustness is this: our ability to obtain nonsingular evolu- 
tion depends importantly on the structure of these equations of motion. Signs are crucial, 
coefficients are crucial, and so is the inclusion of Ramond-Ramond fields. In other words, 
our nonsingular results are highly specific to the field couplings arising from the low-energy 
approximation to string theory. Other "S-branes" arising from "string-motivated" actions 
will probably not possess similarly nonsingular behavior. The positive type of robustness we 
refer to is also a desirable property. What we see manifestly is that the precise form of the 
potentials {V{T),f{T)} is not important, apart from the large-|r| behavior which had been 
derived elsewhere. The most obvious refinement of our work here will be to attack the problem 
of relaxing the requirement of zero NS-NS 5-field. Allowing B^"^^ to be turned on will allow 
us to break ISO{p+l) on the worldvolume and allow inhomogeneous tachyonic modes — the 
importance of which is discussed in refs. |17, 33| — and also to turn on more components of 
R-R fields. Inhomogeneities would of course have to be included in initial conditions, because 
homogeneous on-shell tachyons do not couple to non-homogeneous tachyons |32]. We have 
postponed the non-homogeneous problem to the future mainly because it is messy; our work 
reported here should be considered a step in a larger program. 

The other important approximation we made in our work was in section 3^ , where we had 
to smear the SD-brane sources in the transverse space to facilitate integrating the equations 
numerically by turning them into ODE's. This limits our ability to fully probe the properties 
of the system in which we are interested. Here we would also like to record another physical 
consequence of this ansatz. Namely, this restriction has notable, negative, consequences for 
our ability to track whether black holes form as intermediate states during the time evolution 
of our coupled system including full backreaction. The issue of black holes was raised in 
the discussion section of ref. |^. The essential point is that a black hole intermediate state 
may arise as an alternative to SD-brane formation and decay, at least with the half-advanced, 
half-retarded propagator. The fine-tuned nature of the initial conditions producing SD-branes 
highlights a reason why integrating partial differential equations of motion (including depen- 
dence on transverse coordinates) may be particularly difficult numerically. Or the obstruction 
to finding the full solution may yet turn out to be negotiable. It will also be interesting to 
think further about particle/string production. 

Let us end with some somewhat speculative remarks. Typically in the limit ggN 
the open and closed strings decouple. This is true in our effective lowest-modes analysis here, 
but also explicit in other worldsheet-inspired approaches. There should also exist a limit (in 
time) to be taken where only the open strings survive. From the worldsheet definition of a 
SD-brane, it is suggestive that the open string degrees of freedom would combine to form a 
Euclidean conformal field theory in p+1 dimensions. The same appears true when considering 
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the effective action of massless open string degrees of freedom on an unstable D-brane |34 
In both approaches, however, it is not clear what the role of the tachyon could be. Physically, 
it is the source of a process by which energy is siphoned out of the open string sector and 
pumped into the closed string sector. So, in a sense, the decay of a D-brane through tachyon 
condensation corresponds to the decrease of a c-function-like quantity on the gauge theory 
side. Then, we can entertain the idea that time-evolution on the gravity side should really be 
regarded as a renormalization group (RG) flow on the gauge theory side. Prom this viewpoint, 
formation and decay of a SD-brane would be a process corresponding to first an inverse RG 
flow (integrating in degrees of freedom) followed by regular RG flow (integrating out degrees 
of freedom). This might be related to the study of open string tachyon condensation using 



RG flow in the worldsheet theory |4C]. 
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'Similar ideas were explored in the context of the dS/CFT correspondence (see, for example, refs. 
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A. Regular KMP SD-brane solutions 



Among the supergravity solutions found by there exist some that are regular on the 
horizon located at t = to. This is only realized for the following values of the parameters, 

k = 2, H=^^, G = ki = 0. (A.l) 
7 — p 

The corresponding metric is 



PJ 



(A.2) 



Because these solutions are anisotropic in the worldvolume directions, it is not clear that they 
are physically relevant. We will nevertheless study some interesting properties not considered 
in ref. (5|. For example, the region t = +uj does not correspond to an horizon as suggested 
by the fact that none of the metric components either vanish or diverge there. For p odd the 
solutions are time-reversal symmetric so the region t = —u is also not an horizon. 

A.l The region close to the origin 

For the anisotropic solutions there is no curvature singularity at t = lo. It is therefore 
interesting to consider the behavior of the metric components and curvature invariants close 
to the potentially problematic region around the origin, t = 0. We evaluated the curvature 
invariants: IZ, TZ^yTZ^^ , and Tt^ivpx'Tl^^'^^ . They identically vanish for i = 0. The metric 
tensor there appears suspicious (for example, the component gu diverges) but we find 

P+i 

hm ds^ = -dT^ + T^dHl_p + ^ {dx') \ (A.3) 
^ 1=1 



where r = lo"^ /t. The expression ( A.3 ) is simply flat space with part of it written in Milne 
coordinates. Not surprisingly, we find 

lim $(t) = 0, lim C(t) = 0, (A.4) 
which implies that all stress-energy components vanish in the the region close to the origin. 
A.2 Horizon physics 

We demonstrate that, contrary to previous claims, many of the anisotropic solutions are 
actually regular in the full range: — cxd < t < -|-oo. The anisotropic solutions were already 
shown to be non-singular at t = a; and t = 0. We now investigate the region t = —oj further. 
Let us introduce the change of coordinates 

T= l + (-) t (A.5) 
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in order to make comparison with the results of ref. easier. For p even, T = corresponds 
to t = —io while for p odd we have T = —2io when t = —m. A comment in ref. |^ is that T = 
corresponds to a (non-naked) curvature singularity. Actually, this is not always the case! For 
example, we considered the case p = 1 and computed the associated curvature invariants at 
all times. Figure |^ shows the evolution of the Ricci scalar for the solution with co = 1 and 
9 = 7r/4. Clearly, the p = 1 solution is symmetric under time-reversal and therefore has no 
curvature singularity. It also does not possess any horizon, a feature common to the regular 
solutions found in this paper. The qualitative behavior of all curvature invariants is similar 
to the Ricci scalar and is quite generic, i.e., it is unchanged for all odd values of p, 9 and uj. 
For p = 1 we obtain 

32^/3 19 cos^ 61 - 26 cos2 (9 + 7 , 
lim7^ = ^ E . (A. 6) 

This is finite except for = in which case the Ricci scalar diverges like R ~ — iSf" . We 
also found expressions for two other curvature invariants (for p = 1), 

_ 923/4-272cos6 6l + 14 + 255cos''6l + 101cos'^6l-98cos2(9 
t^±^^ 2u;^ sm^ 9 

322/3-ll2cos6l + 28 + 892cos8 6l- 1840cos6 6l + 1032cosS6l 
lim IZ^ypalZ^ - — ri^-- . (A.8) 

We found expressions with similar qualitative behavior for other values of p odd. 

For p even the solutions are not time-reversal symmetric. As pointed out previously, 
curvature invariants are finite at t = a; but there is a curvature singularity at t = —a;. These 
are the spacelike curvature singularities (protected by an horizon at t = 0) described in ref. |p. 

A. 3 Unstable brane probe analysis 



As mentioned in section 2.2, the motivation behind considering an unstable brane probe in a 
background with singularity problems is to ask if the singular background could actually be 
built. 

The calculations and results of ref. were summarized in section p.2| . In this appendix 
we generalize this calculation by probing the d = 10 anisotropic backgrounds presented above, 
since these are the only ones which are either non-singular or have singularities (at t = —to) 
shielded by a horizon (t = 0). We felt this generalization to be necessary because ref. did 
not, for example, address the issue as to how the inclusion of the dilaton might affect the 
brane probe calculation. The unstable brane action is the obvious generalization eq. ( ^.61) of 
the case studied in ref. [0. 

We investigate whether or not an unstable brane probe is a well-defined object in the 
vicinity of the region t = uj. In Einstein frame, the energy density for the probe propagating 
in the anisotropic backgrounds is 

P,ro,e = ^m)^, (A.9) 
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Figure 7: This figure illustrates the Ricci scalar for the p = 1 anisotropic SD-branc solution with 
Lu = 1 and 6 — tt/A. The other curvature invariants behave similarly. 



while the pressure corresponds to 



Pprobe = ^^/($)y(r)AV2. (A.IO) 

gs 

The dilaton function /($) was picked up during the transformation from the string frame 
to the Einstein frame. It plays no role in the upcoming analysis because the dilaton is well- 
behaved, 

lim /(^>) = const. (A.ll) 



As we saw in section 2.2, whenever the probe analysis goes wrong, it signals a pathology 
for the gravitational background. As we saw, there are at least two ways the probe analysis 
can go wrong: it may induce an infinite energy or pressure density {pprobe^Pprobe ioo), or, 
there might not exist any reasonable solutions for T(t). 

For t LO, the dominant contribution to the equation of motion for the tachyon is 

9 2 

A2-A-F-(t-u;)A = 0. (A.12) 

This is solved for 

[t-ujyi^ - g 

where 5 is a constant of integration. The solution A = (5 7^ 0) clearly corresponds to the 
brane probe inducing a curvature singularity on the horizon. The only physical solution is the 
one for which g = which corresponds to A = 1. For the anisotropic backgrounds considered 
here the metric component gu neither vanishes nor blows up at t = to. The solution A = 1 
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therefore corresponds to a tachyon field for which the time-derivative vanishes (T = 0) at 
t = uj. It therefore appears that there are solutions for the probe evolution that avoids the 
pathologies described earlier. This is no surprise since for these anisotropic backgrounds the 
region t = uj is not an horizon in the technical sense of the term. We repeated the calculation 
around the regions t = and t = —uj. We find that for both p odd and even the unstable 
brane probe is not well-behaved at t = 0, i.e., it induces a curvature singularity there. 



B. Tachyon in flat space 

We consider solutions to the equation of motion for an open string tachyon when the massless 
closed string modes are decoupled. The relevant equation of motion is 

B.l General solution 

Eq. ( |B.lD is a second order differential equation with a general solution of the form 

where a and b are constants of integration. To integrate this equation one needs the function 
V{t) which would imply that we already know the solution T{t). Open string field theory has 
taught us that for t = tc large we have 

lim V{t) « 1, (B.3) 

t — >tc 

which implies that 

T{t) ~ / dt{l + 2aV(t)) + b. (B.4) 

Jtc 

Therefore at large time the tachyon behaves like 

T{t) = t + {b-tc) + 2c? [ dt V^{t). (B.5) 

Jtc 

Using the string field theory result 

lim V(T) = e"^/^, (B.6) 

T-^+oo 

and taking T{t) ~ t leads to the large time formula 

T{t) ~ t - a'^V2e-^\ (B.7) 
where we have fixed the integration constants by imposing 

a^V2e-^*^ + b-tc = 0. (B.8) 
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B.2 Particular solutions 

We consider solutions to the equation of motion (|B.1D with the potential 



The tachyon equation of motion becomes 

r+ -^(1 - T2)tanh (^r/y2) , (B.IO) 
which has a solution of the form 

T{t) = -^/2arcsinh [cie*/^ - cae"*/^] ^ , (B.ll) 

where c\ and C2 are constants of integration. We usually specify boundary conditions at t = 0, 

r(0) = -72 arcsinh |^ ^/^(ci^- C2) ^ ^^^2) 

r(0) = -V2 "-^^ ^. (B.13) 

The family of solutions characterized by T(0) = (ci = C2) corresponds to all possible 
tachyon velocities at t = 0: T(0) = —\plc\. An other class of solutions are those for which 
T(0) = (ci = — C2). Those correspond to allowing the tachyon to begin its evolution with 
T(0) = -\/2 arcsinh ^/2ci. 

The solution we presented are referred to as tachyon matter. The stress-energy compo- 
nents (which are independent of the number of dimensions in the theory) correspond to a 
conserved energy density (p ~ y(T)/\/A = constant) and a pressure ~ —ViT^^fK) that 
vanishes as t — > +00 [llSll. 
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